Nodal Domain Statistics for Quantum Maps, Percolation and SLE 
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We develop a percolation model for nodal domains in the eigenvectors of quantum chaotic torus 
maps. Our model follows directly from the assumption that the quantum maps are described 
by random matrix theory. Its accuracy in predicting statistical properties of the nodal domains is 
demonstrated by numerical computations for perturbed cat maps and supports the use of percolation 
theory to describe the wave functions of general hamiltonian systems, where the validity of the 
underlying assumptions is much less clear. We also demonstrate that the nodal domains of the 
perturbed cat maps obey the Cardy crossing formula and find evidence that the boundaries of 
the nodal domains are described by SLE with k close to the expected value of 6, suggesting that 
quantum chaotic wave functions may exhibit conformal invariance in the semiclassical limit. 
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One of the central problems in the field of quantum 
chaos is to understand the morphology of quantum eigen- 
functions in classically chaotic systems. In time-reversal- 
symmetric systems one can always find a basis in which 
these eigenfunctions are real. They can thus be divided 
into nodal domains - connected regions of the same sign, 
separated by nodal lines on which the eigenfunctions van- 
ish. The statistical properties of these nodal domains 
then constitute a natural way to characterize the mor- 
phology of the eigenfunctions. 

Nodal domain statistics were studied for separable bil- 
liards in 0, where it was shown that if v n is the num- 
ber of nodal domains in the nth energy eigenstate then 
Xn = v n jn has a limiting distribution as n — > oo with a 
square-root singularity at a system-dependent maximum 

Value Xmax < 1- 

In chaotic systems the eigenfunctions may be modeled 
statistically, far from boundaries and turning points, by 
random superpositions of plane waves 0: 



u(x) 



cos(kx cos 0j + ky sin 9j + <fij) (1) 



where 9j and 4>j are random phases. This is known 
as the random wave model. Since plane waves are so- 
lutions of the Schrodinger equation for a free particle, 
V 2 ?/> = —k 2 ip, the maxima of any superposition are pos- 
itive and the minima are negative. Hence the nodal do- 
mains correspond to groups of either maxima or minima. 
A given pair of adjacent maxima (minima) lie in the same 
nodal domain if the saddle point between them is pos- 
itive (negative). The density of saddles in the random 
wave model is asymptotically twice the density maxima 
or minima. This would be exactly the case, for exam- 
ple, if the maxima and minima lay on alternate sites of a 
square lattice and the saddles on the corresponding dual 
lattice, e.g. midway between diagonally adjacent maxima 
(or minima) (although it is important to note that in the 
random wave model typical realizations of the wave func- 
tion are in fact highly irregular). The saddles may be 



thought of as lying at the midpoints of bonds of the dual 
lattice connecting the maxima, for example. If the sad- 
dle height is positive, then the corresponding maxima are 
connected and the bond may be thought of as 'open'; if it 
is negative, the maxima are not directly connected, and 
the bond may be thought of as 'closed'. This was the 
basis of the very interesting suggestion put forward by 
Bogomolny and Schmit Q that statistical properties of 
nodal domains in the random wave model, and hence in 
quantum chaotic eigenfunctions, correspond to those in 
critical percolation. Specifically, Bogomolny and Schmit 
assumed that the heights of the saddles are uncorrelated 
and have equal probability of being positive or negative, 
and proposed bond percolation on a square lattice as a 
model for nodal domain statistics. This implies that Xn 
is Gaussian distributed as n varies in the semiclassical 
limit. Moreover, it leads to the conclusion that the scal- 
ing exponents associated with critical percolation also 
characterize properties of the nodal domains in quantum 
chaotic eigenfunctions, for example their area distribu- 
tion and fractal dimension. 

The predictions of the percolation model are consistent 
with numerical computations [2| and experimental mea- 
surements |4j of the fluctuation statistics for the nodal 
domains of quantum billiards, but the data do not pro- 
vide conclusive verification. This is important, because 
the model has been the subject of considerable debate. 
Foltin has shown that the heights of the saddles in the 
random wave model exhibit long range correlations 0|, 
contradicting one of the key assumptions of the perco- 
lation model. Bogomolny has argued that oscillations in 
the two-point correlation function are sufficient to ensure 
the applicability of the Harris criterion and so guar- 
antee that the scaling exponents are unaffected Q, but 
the issue awaits a more systematic investigation. More- 
over, Foltin, Gnutzmann and Smilansky have devised a 
particular statistical measure for which the percolation 
model fails The range of validity of the model and 
the precise assumptions upon which it relies thus remain 
to be determined. 
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Our purpose here is to establish a percolation model 
for quantum torus maps. These are some of the most im- 
portant examples of quantum chaotic systems, because 
one can find maps that are fully chaotic and quantum 
mechanically they are finite dimensional and so easily 
tractable. We will show that for these systems there is 
a critical percolation model that follows directly from 
the Bohigas-Giannoni-Schmit (BGS) conjecture, which 
asserts that local quantum fluctuation statistics in clas- 
sically chaotic systems are modeled by Random Matrix 
Theory This model corresponds to site percolation on 
a triangular lattice, which falls into the same universality 
class as bond percolation on a square lattice and so has 
the same critical exponents. The advantages of investi- 
gating the percolation model for maps are, first, that the 
assumptions underlying it are very much more straight- 
forward - one only has to assume the BGS conjecture, 
and there are no problems analogous to those relating 
to the slow decay of correlations in billiard eigenfunc- 
tions - and, second, that one can perform more exten- 
sive and controlled numerical computations, leading to 
significantly more precise tests of the predictions. 

We find that the percolation model for maps is ex- 
tremely accurate in that the critical scaling exponents 
associated with the nodal domains are very close to those 
predicted by percolation theory. Moreover, the agree- 
ment goes beyond scaling laws: the nodal domains of the 
quantum maps we study also obey the Cardy crossing for- 
mula and its generalizations. We also verify that Cardy's 
formula is satisfied within the random wave model. This 
suggests that both linear superpositions of random waves 
and quantum chaotic eigenfunctions may exhibit confor- 
mal invariance in the semiclassical limit. Finally, the 
link between processes governed by Stochastic Loewner 
Evolution (SLE) and statistical models has recently been 
the focus of considerable attention. Critical percolation 
is believed to relate to SLE with diffusion constant n = 6. 
For percolation on a triangular lattice this has been es- 
tablished rigorously ■ O n the basis of the percolation 
model one would expect nodal lines to behave like pro- 
cesses governed by SLE 6 . We find evidence that this is 
the case for the boundaries of the nodal domains in the 
case of quantum maps. 

The systems we study correspond to chaotic symplctic 
maps acting on the unit 2L-dimensional torus, which is 
viewed as their phase space. Such maps can be quan- 
tized using an approach introduced by Hannay and Berry 
(llj . The Hilbert space has finite dimension N L , where 
N plays the role of the inverse of Planck's constant. 
Quantum maps correspond to unitary matrices U act- 
ing on wave functions in this Hilbert space so as to gen- 
erate their (discrete) time evolution. In the position 
representation these wave functions take values on an 
L-dimcnsional lattice. For example, when L = 1 the 
wave functions take values (ci, C2, ■ ■ • , cjsr) at positions 
q = Q/N, < Q < N; and when L = 2 they take 



values (ci, ca, . . • , Cya) at positions on the square lattice 
q = (Q 1 /N,Q 2 /N), < Q U Q 2 < N. We shall be con- 
cerned with the quantum map eigenvectors. If the map 
is time- reversal symmetric (and so U is symmetric), the 
components of the eigenvectors are real. For a given 
eigenvector, we can thus split the quantum lattice into 
regions such that the components associated with neigh- 
bouring sites have the same sign. These regions then 
correspond to nodal domains. 

When L = 1 this can be done straightforwardly: if sites 
lying next to each other on the 1-dimensional lattice have 
eigenvector components Cj with the same sign then they 
constitute part of the same nodal domain. When L = 2 
the situation requires more careful consideration, because 
one needs a convention for deciding whether lattice points 
that are diagonal neighbours and have eigenvector com- 
ponents with the same sign lie in the same nodal do- 
main or not. Consider, for example, when the eigenvec- 
tor components associated with a group of four lattice 
points which form a square have signs in a checkerboard 
arrangement, e.g. on the top row + -, and underneath - 
+. Are the pluses automatically part of the same nodal 
domain, or the minuses? We take as our convention that 
lattice points are connected along diagonals running from 
the top left to the bottom right; so in the example just 
given it is the pluses that are connected. This takes us 
from the original square lattice to the triangular lattice. 
Nodal domains then correspond to regions on this tri- 
angular lattice in which connected points have the same 
sign. Our convention is, of course, one of many possibili- 
ties. However, we note that it is necessary to incorporate 
diagonal neighbours for the definition of nodal domains 
to be consistent with that in billiards, and that all of 
the conventions we have tested which do this lead to the 
same results. 

In order to develop a statistical model for the nodal 
domains we now need to introduce a statistical ansatz 
for the signs. According to the BGS conjecture, for 
generic, classically chaotic, time-reversal-symmetric sys- 
tems statistical properties of the matrix U should coin- 
cide with those of random matrices taken from the Cir- 
cular Orthogonal Ensemble of Random Matrix Theory. 
The joint probability density for the eigenvector compo- 
nents (c% , C2, . . . , c/v-l ) is then 

P( Cl ,c 2 ,...,c NL ) = S ^$>f-lj (2) 

that is, the eigenvectors are uniformly distributed on the 
unit hypersphere. Crucially for us, it follows immediately 
that the sign of a given component is equally likely to be 
positive or negative and that these signs are independent 
of each other at different sites, i.e. they are uncorrelated. 

When L = 1 this model was explored in [J^. When 
L = 2 it corresponds directly to site percolation on a 
triangular lattice, which falls into the same universality 
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class as the Bogomolny-Schmit model. This means that 
the critical exponents associated with the nodal domain 
statistics will be the same. We note that in our case these 
have been established rigorously for percolation |icj . 

We now test the percolation model for a particular 
family of quantum torus maps. In essence, we are seeing 
whether this family is described sufficiently accurately 
by RMT for the model to apply. Linear maps are not 
sufficient for our purpose: because of non-generic arith- 
metical symmetries they are not described by RMT ^jj . 
Instead, we take a linear map composed with a nonlinear 
perturbation. Specifically, we use M — p o Ao p with 
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(3) 

and p a nonlinear periodic shear in momentum: p\ — » 
pi + fi cos(27rgi), p 2 -» pi + |j cos(27rg 2 ). The map M 
is time-reversal-symmetric and, for sufficiently small val- 
ues of the perturbation parameters, completely chaotic. 
The corresponding quantum map, a unitary matrix of 
dimension N 2 can be written down easily using the pre- 
scriptions in [13 and [3 (the formula is long and so we 
do not give it here, see [l5| for details). We now compare 
statistical properties of the nodal domains of this map 
with those of percolation clusters. 

Consider first the distribution of the number n of nodal 
domains. For percolation on iV 2 sites this should be a 
Gaussian with mean n c N 2 + b + o(l) and variance cN 2 , 
where Monte Carlo simulations give n c = 0.0176 . . ., b — 
0.878 ... and c = 0.0309 . . . 16]. For the quantum map 
we find a Gaussian distribution with a mean and variance 
consistent with the percolation formulae. Specifically, 
when k% — 0.04 and k 2 = 0.01, a best fit gives n c = 
0.0176, b = 0.902 and c = 0.0297. The data are shown 
in Figure ^ For the distribution of areas a of the nodal 
domains, the percolation model predicts a scaling law 
a~ T , with t = 187/91. A log-log plot of the data for 
the eigenvectors of the quantum map is shown in Figure 
|5J The percolation model also implies that the nodal 
domains should have a fractal dimension D = 91/48 = 
1.89.... Data for the quantum map, obtained using a 
box-counting algorithm and shown as a log-log plot in 
Figure |31 are consistent with this in that the best fitting 
straight line has a gradient 1.8774. 

One of the key features of critical percolation is that it 
has a conformally invariant limit. This underlies the use 
of conformal field theory, in deriving the Cardy crossing 
formula for example, and the link with SLE. Given the 
success of the percolation model in describing scaling ex- 
ponents associated with their nodal domains, it is natural 
to ask whether random waves and quantum eigenfunc- 
tions are also conformally invariant in the semiclassical 
limit. In percolation, the Cardy crossing formula gives 
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FIG. 1: The mean (red circles) and variance (blue squares) 
of the number of nodal domains for the quantum map with 
fei = 0.04 and ki = 0.01. The linear fit for the mean (solid 
line) gives n c — 0.0176 and b = 0.902, while the fit for the 
variance (dashed line) gives c = 0.0297. 




FIG. 2: Frequency of nodal domains as a function of area. 
The solid (red) line shows the data and the dashed (blue) line 
the theoretical power law with exponent r = 187/91. Here 
N = 61, fci = 0.04 and k 2 = 0.01. 



the probability of there being a cluster spanning the sys- 
tem between fixed sections of the boundary [17]. For 
the eigenvectors of the quantum map, we can determine 
the probability that a nodal domain spans the lattice, 
connecting a fixed section of length x on the left to any 
region on the right-hand side. We also test the nodal do- 
mains of realizations of the random wave model using the 
same geometry. The results are shown in Figure where 
they are compared to the crossing formula. To explore 
the connection with SLE we use an idea introduced in 
(lif to test conformal invariance in 2-D turbulence. The 
SLE K process generates a trace from a stochastic differen- 
tial equation dependent on a driving function £ (t) , which 
is proportional to a Brownian motion. The trace is ob- 
tained for each eigenvector by following a nodal line while 
keeping the positive points to the right. Upon hitting the 
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FIG. 3: A box count for the largest nodal domain of an eigen- 
vector of the map with TV = 61, fci = 0.01 and k 2 = 0.02. The 
linear fit corresponds to a fractal dimension of 1.8774. 




FIG. 4: Numerically computed crossing probabilities for map 
eigenvectors when TV = 61, fci = 0.01 and k 2 = 0.02 (red 
diamonds) and 2000 realizations of the random wave model 
(blue squares), compared with Cardy's formula (solid line). 

left boundary the trace continues along the boundary, in 
such a way that it can always be connected to the right 
boundary without crossing itself, until it reaches another 
nodal line. The process is stopped when the trace reaches 
any of the other boundaries. The driving is computed by 
approximating the candidate trace by a union of straight 
line sections. The exact solution of the Loewner evolu- 
tion for a straight line growing from the boundary [19j is 
then used to find maps which 'swallow' each straight line 
section in turn. This gives a sequence of times and driv- 
ing parameters. The results are shown in Figure [SJ The 
Gaussian distribution (inset) is consistent with Brownian 
motion. The best fit to the diffusion constant is k ks 5.3. 
This method is self consistent because we obtain a value 
of kappa close to 6. SLE6 has the locality property that 
the trace does not feel the effect of the boundary until 
it touches it. This allows us to use the half- plane SLE 
process even in this restricted geometry. 

The fact that the crossing formula applies and the link 
with SLE holds is evidence of conformal invariance. It 
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FIG. 5: Variance of as a function of t. The inset shows the 
distribution of the rescaled driving function for t = 0.3 (red 
circles), t — 0.4 (blue squares) and t — 0.45 (green diamonds) 
with the expected gaussian shown as a solid line. Here N = 
61, fci = 0.02 and k 2 = 0.01. 



is natural to conjecture that this will extend to generic 
quantum chaotic eigenfunctions in the semiclassical limit. 
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